Quantum Depletion of a Soliton Condensate 
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Abstract 
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q ■ We present rigorous results on the diagonalization of Bogoliubov Hamilto- 

nian for a soliton condensate. Using the complete and orthonomalized set 
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l/~j ■ of eigenfunction for the Bogoliubov de Gennes equations, we calculate ex- 

actly the quantum depletion of the condensate and show that two degenerate 
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vq ■ zero-modes, which originate from a U(l) guage- and a translational-symmetry 
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breaking of the system, induce the quantum diffusion and transverse instabil- 
ity of the soliton condensate. 
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The study of matter waves and elementary excitations has received much attention be- 
cause of the remarkable experimental realization of Bose- Einstein condensation in trapped, 
weakly interacting atomic gases [1-3]. In the past few years, several paths have been ex- 
plored in studying elementary excitations in Bose-Einstein condensates (BECs). The most 
followed method is to use the Gross-Pitaevskii (GP) equation which is suitable to describe a 
zero-temperature BEC One of the deficiency of this approach, in addition to the constraint 
of zero temperature, is the neglect of the quantum fluctuations of condensate, which is an 
important aspect of elementary excitation [4]. The eigen-modes of the elementary excita- 
tions thus obtained so far [1,2] are not complete and their orthonormalities have never been 
proved. 

A different approach on elementary excitations is to use Bogoliubov theory [5] that was 
originally formulated for homogeneous Bose systems but is also valid for inhomogeneous ones. 
In this approach one makes a canonical transformation for boson operators to diagonalize 
the quantum Hamiltonian of system and hence the quantum fluctuations of condensate are 
taken into account [6]. An added advantage of the Bogoliubov theory is that it can be 
easily generalized to the case of finite temperature (with, e. g., the Hatree-Fock and Popov 
approximations) [1,2], therefore can be compared directly with experiments carried out at 
non-zero temperature environment. 

In this Letter we apply the Bogoliubov theory to the investigation of the elementary 
excitations generated from a quasi one-dimensional (ID) soliton condensate. It is known that 
a quasi ID BEC can form in the presence of a transverse confinement, which introduces a 
cutoff for long wavelength fluctuations in the transverse directions and hence an off-diagonal- 
long-range order can establish, as have been reported in the case of low-dimensions BEC 
systems [7] and in the case of matter- wave bright solitons in quasi- ID systems [8,9] and 
the corresponding flourished theoretical activities [10-13]. The main contribution of our 
present work is the rigorous diagonalization of quantum Bogoliubov Hamiltonian and the 
exact calculation of quantum depletion of the soliton condensate and related conclusions on 
soliton stability. 



We consider a quasi ID attractive, weakly interacting Bose gas condensed in a trap with 
tight transverse confinement and negligible trapping potential in the axial direction. The 
grand canonical Hamiltonian of the system can be written as the following dimensionless 
from 

d 2 



K = H-fiN = fdz&(z,t) 



H-W(z,t)V(z,t) 



*(*,*), (1) 



dz 2 

where ^f(z, t) is the field operator satisfying [&(z, t), W(z' , t)] = S(z-z') (with other commu- 
tators being zero). The chemical potential \x assures the conservation of the average particle 
number N = J dz^^(z,t)^>(z,t). In the expression (1) the coordinate, time, field operator, 
and energy (including chemical potential) are scaled by / = h 2 /(m\Uo\), 2ml 2 /h, l~ 1 ^ 2 , and 
h /(2ml 2 ), respectively. U$(< 0) is the interatomic interaction constant and m is the atomic 
mass. The model (1) has been used in Ref. [12] to study the quantum phase transition 
from a homogeneous ground state to a bright-soliton ground state when the strength of the 
attractive interaction is increased. 

Letting ^(z,t) = ip g (z) + ip(z,t) and assuming \\ft\ << ip g we obtain the Bogoliubov 
form (quadratic for ip and ft) of the grand canonical Hamiltonian K = K g + K 2 with K g = 
J dzip g (-d 2 /dz 2 -fi- ip 2 )ilj g and K 2 = J dz [^(-d 2 /dz 2 - // - 4^)^ - ip 2 g (^ + ft ft)] , 
where the ground-state wave function satisfies the equation (—d 2 /dz 2 — \x — 2ip 2 )if) g = 0. For 
small \Uq\ the ground state is a homogeneous condensate and hence the Hamiltonian can 
be easily diagonalized and the elementary excitation exhibits a gapless phonon spectrum. 
For large \Uq\, however, this homogeneous condensate is dynamically unstable [12] and the 
system undergoes a quantum phase transition into the soliton ground state 

i, g {z) = (iV /2)sech[(iV /2)(^ - z )} (2) 

with chemical potential \x = —Nq/4, where Zq is an arbitrary constant, indicating that the 
quantum phase transition results in a translational symmetry breaking of the system. 

Our first objective is to search for a rigorous diagonalization of the Hamiltonian of the 
soliton condensate (2). Letting (N /2)(z — z ) — > z and t/j = (iVo/2) 1 / 2 we get 



K = K g + (iV 2 /4) f dz [$1$ - sech 2 z(00 + f f )] , (3) 

where L = —d 2 /dz 2 — 4sech z + 1 and <p satisfies [0(z, £), $(z', t)] = S(z — z') (with other 
commutators being zero). To diagonalize K we make the canonical transformation 



J2 [uj(z)cj + v*{z)c) 



+ J2 [M z )ck + v* k {z)c\ 



(4) 



where dj and c k are usual boson operators. The first and second terms on the right-hand- 
side (RHS) of (4) are respectively the contributions from the discrete (j) and continuum (k) 
spectra of the excitations generated from the soliton condensate. The key for diagonalizing 
K is to find a complete set of the eigen-functions {u q (z),v q (z);q = k,j} with which the 
expansion of <fi can be made. The discrete modes are in fact the eigen-modes with a zero 
eigenvalue, as shown below. 

Substituting (4) into (3) and assuming that u q and v q fulfill the following eigenequations: 

L 2 L l( t> q {z) = E 2 <f> q (z), L x L^ q (z) = E 2 ^ q (z), (5) 

where ip q (z) = u q (z) +v q (z), 4> q (z) = u q (z) —v q (z) and Lj = d 2 /dz 2 + (25 j\ + 6^- 2 )sech 2 z — 1 
(j = 1, 2), by a detailed calculation we obtain a diagonalized K = K g + (Nq/4)K 2 o, with 

(6) 



K 20 = --J2 [ dzE k \v k (z)\ 2 -P 2 + Q 2 2 + J2 E kclc k , 
6 '-^ J fc^O 



fe^O' 



where we have introduced the operators [6] Pj = (6j + cj)/\/2 and Qj 



i(cj - c])/V2. 



Equation (5) is equivalent to the Bogoliubov-de Gennes (BdG) eigen-value problem: 



Lu q (z) — 2sech 2 ,2t> g (,2) = E q u q (z), 



Lv q (z) — 2sech zu q (z) 



■E q v q {z). 



(7) 
(8) 



The solutions of the Eq. (5) has been known in classical soliton perturbation theory [14,15]. 
Thus we can use the result obtained in Refs. [14,15] to get the following solutions of the 
BdG equations (7) and (8): 

1 



u k (z 



Vk(Z 



/2tt(/c 2 + 1; 

1 
^7r(fc 2 + 1) 



k + 2ik tanh z — 2 tanh z exp(ikz) 
sech zexp(ikz) 



(9) 
(10) 



with 

E k = k 2 + l (11) 

for k 7^ (continuous spectrum). Equations (7) and (8) admit also the following zero-mode 
solutions (j = 1, 2): 

Uj(z) = Sji sechz (2 — ztanhz)/2 + Sj2 sechz [z + tanhz)/2, (12) 

Vj(z) = —Sji (z/2) sechz tanhz + Sj2 sechz (tanhz — z)j1. (13) 

These zero-modes belong to the discrete spectrum of the system with a zero eigenvalue, 
which is two-fold degenerate. Physically, such degenerate zero-modes originate from a U(l) 
guage- as well as a translational-symmetry breaking of the system and thus open a gap in 
the excitation spectrum (see Eq. (11) ). This is very different from the result obtained for 
homogeneous condensates. It can be shown that the eigenfunction set obtained above form 
a complete function set and they are also orthogonal [16]. 

With the above results we find that K g = 2Nq/3 and J^k^o^k f dz\v(z, k)\ 2 
1)/(6tt) and obtain diagonalized (dimensional) Bogoliubov quantum Hamiltonian 



[71 - 



24 ft 2 ° 8 h 2 ° 



6 D7r fc^O 



(14) 



Note that each term on the RHS of (14) has clear physical meaning. The first term describes 
the ground state energy of the system, which for large N agrees with the result of Bethe 
ansatz solution [17]. The first term in the square bracket comes from the zero- modes whereas 
the second term is the contribution of the depletion of the condensate. The terms —P 2 + Q\ 
are originated from the quantum fluctuations of the condensate and the last term represents 
the contribution from phonons. Because the continuous spectrum of the excitations E k 
for the soliton condensate opens a gap, thus the phonons are massive due the breaking of 
translational symmetry of the system. This is very different from the case of homogeneous 
condensates where a phonon is non- massive. 

Our next objective is to investigate the quantum dynamics of the soliton condensate. 
For the zero-modes we obtain the equations of motion 



dP 1 /dr = 0, dQ 1 /dr = -2P U 

dP 2 /dr = -2Q 2 , dQ 2 /dr = 0, 



(15) 
(16) 



where r = (iVg/4)£. The Hermitian operators Pj and Qj satisfy the commutation relation 
[Qj, Pj'] = iSjj', and they are associated with the collective motion of the solition condensate 
[6]. The exact solutions of the above equations read Piij) = A(0), Qiij) = Qi(0)— 2P 1 (0)r, 
A(t) = -£2(0) — 2Q 2 (0)t, and Q 2 {r) = Q 2 (0). Using the zero-mode solutions and the 



relations between dj, dj and Pj, Qj we obtain 



$(z,t) ps A(z,t) sech 



z — 



Wn 



--Pi.(t) 



cxp 



"No 



Qi(r)+zQ 2 (r) 



(17) 



with A(z,t) ps (Nq/2) 1 + (1 — z tanh z) Pi(r) / \^1% . Equation (17) clearly shows quan- 
tum fluctuations in the amplitude, position and phase of the soliton condensate due to the 
existence of the zero-modes. This is again very different from the case of a homogeneous 
condensate. 

The above results can be used to calculate the quantum depletion of the soliton conden- 
sate [18]. Let \N , Ni, N 2 , N exc > denote the state with A^ atoms in the condensate. JVj 
and N 2 are number of atoms occupying the first and the second zero-modes, and N exc are 
the number of atoms being in the continuous modes. Assume that the initial (t = 0) state 
of the system is a Bogoliubov (quasiparticle) vacuum \N , 0,0,0 >. We consider the time 
evolution of the particle numbers occupying in different modes. Noting that the particle 
number density operator of the system is given by n(z,t) = l~ 1 ^(z,t)^f(z,t) and using the 
expression of the field operator, we get the particle-number density for the condensate mode 
n (z)=< rio(z) >= l~ 1 (N /2) 2 sech z, which gives rise naturally the total particle number 
Nq in the condensate. Here < • • • > denotes the average over the Bogoliubov vacuum. It is 
easy to show that the particle-number densities in the first and the second zero-modes are 
given by 



nj(z,r) = ^ [v](z) + [uj(z) + (-iy V j(z)] 2 r 2 } , (j = 1,2). 



(18) 



n exc (z, t) = -^J2 \v k (z) | 2 = -^ (- - —J sechS. (19) 



The atoms in these zero-modes are incoherent ones, representing the quantum depletion of 
the soliton condensate. There are two features for the quantum depletion: (i) The spatial 
distributions of the incoherent atoms in the zero-modes are localized because Uj(z) and Vj(z) 
are localized functions of z; (ii) The distribution densities of the incoherent atoms are time- 
dependent and proportional to r 2 , which implies that the depleted atoms grow algebraically 
as time increases and hence the soliton condensate loses atoms spontaneously even in the 
absence of thermal cloud. This is a manifestation of the quantum diffusion of the soliton 
condensate implied in (17). The exact expressions of the total particle numbers depleted 
in the first and the second zero mode are N\(t)= f^dzni(z,t) =(12 + 7r)/72 + 2r 2 and 
N2(t)=f^ XJ dzn 2 (z,t) =(2/3)(l + r 2 ). The particle-number density in all continuous modes 
is given by 

A-'n^. . ., An /I 1 

Thus, the total particle number depleted in the continuous modes is A^ e;rc (t)=A^ exc (0) = 
J^° QO dzn exc (z,t) = ('7i — 2)/(37r), a small number in comparison with that depleted in the 
zero modes. 

We can also obtain the exact expressions of the average values for Pj and Qj. They are 
given by < A(t)A(t) >=l/2, < p2(r)P 2 (r) >=(l/2)(l + 4r 2 ), < Q 2 (r)Q 2 (r) >=(1/2)(1 + 
4r 2 ) and < Q2(t)Q 2 {t) >=l/2. These expressions reflect again the diffusion property of 
the soliton condensate. It should be noted that the instability of the soliton condensate 
shown above is based on a linear Bogoliubov approach. The inclusion of the nonlinear effect 
between the quasiparticle excitations may modify these results. 

Our third objective is to study the stability of the soliton condensate under a long 
wavelength transverse perturbation. Assuming the confinement in the x and y direc- 
tions is relaxed the operator d 2 /dz 2 in (1) and (3) is replaced by d 2 /dx 2 + d 2 /dy 2 + 
d 2 /dz 2 . To diagonalize the Hamiltonian in this case we use the canonical transformation 
Ek ± Eg u q ,k ± (z) exp(ik ± •r_ L )c g , k± + v* k± (z) exp(-ik ± • rjjc^ , where r± = (x,y) and 
q = (k x , k y , q)—(\a±, q) with q = k (for continuous modes) or j (for discrete modes). The 



eigen equations (5) are replaced by 

(L 2 - klXL, - fci)0 q (z) = ^ q q (z), (20) 

(Li - kl)(L 2 - kl)^(z) = crMz), (21) 

where <r q = E"l. We consider a long wavelength perturbation in the transverse (x and y) 
directions, thus k\ is a small quantity which can be taken as a small parameter to solve 
Eqs. (20) and (21) using a perturbation theory. Taking q = q °) + k\<\)^> + k\<\>^ + • • •, 

■0 q = V'q^ +^l'0q 1 ^ +^i'0q 2 ^ H an d °"q = o~^ + k\a^ +k\a^ H , we get a hierarchy of the 

equations on q '\ ^® and <r q ^ (I = 0, 1, 2, ...). The leading-order solutions, including the zero- 
modes and the continuous modes, are just those obtained given above. We are interested 
in what will happen for the zero-modes (a^> = 0) when the transverse perturbation is 
applied to the system. Thus in the following we consider the dynamics of one zero-mode 
(say (p^(z) = 4>i(z) = Ui(z) — Vi(z) = sechz) [19]. Because the leading-order solutions are 
a set of complete and orthonormalized functions, all high-order solutions can be obtained 
exactly. In the orders I = 1 and 2, for frequency correction we get 0~^——2j^ oo (f)\( K z) = 
-4 and 4 2 ) = 1 + JZo dka[ l) JZo dz [(k 2 + l)^{z)^ k {z) - 2<t> 1 {z)4> k {z)) = 3/2. Here a k 1] = 



— Jir/2k/[(k 2 + l) 2 sinh(7r/c/2)]. Thus under the long wavelength transverse perturbation 
the eigenfrequency corresponding to the zero-mode reads 



E 1M± = ^-4k 2 ± + 3&i/2 = i2k L ^l - 3k 2 j8, (22) 

which is pure imaginary for k\ < 8/3 and hence the time growth rate of the zero-mode is 



Im(i? 1) k x )=2A;_|_w 1 — 3k"j_/8. Another zero-mode displays also similar behavior [16]. Con- 
sequently, if there is no confinement in x and y directions, the soliton condensate will be 
dynamically unstable under the long wavelength (i. e. for small k±) transverse perturbation. 
The origin of this instability comes from the zero-modes. Thus to repress this instability a 
confinement in the x and y directions is necessary. This is why for observing a stable soliton 
condensate in experiment one must use a quasi ID trapping potential. 

We now briefly discuss the experimental aspects of quasi ID condensate. Low dimen- 
sional BECs have been realized experimentally [7] and the dark solitons in the BECs have 

8 



been investigated intensively [20-22] . Bright solitons and bright soliton trains for attractive 
atomic interaction have been observed recently [8,9]. These works will further stimulate in- 
tensive studies on elementary excitations in inhomogeneous Bose systems. Experimentally, 
a quasi-lD BEC, as those suitable for the present work, can be obtained by tightening the 
confinement in two transverse directions so that the energy-level spacings in those directions 
exceed the interaction energy. Consequently, the motion of the atoms in the transverse di- 
rections is frozen out and the system is essentially one dimensional. Note that the results 
presented above can be applied to a quasi-lD Bose gas with attractive atomic interaction. 
By assuming the 3D field operator with the form $(r, t)=§o(r±)ty(z,t), where $o( r ±) is 
the single-particle ground state wave function in the transverse directions, one can obtain a 
quasi-lD Hamiltonian with the same form of (1) but the coordinate z should be scaled by 
/ = /*7(m/ |^|) with U' Q = £/ JdrJ$o(r ± )| 4 . 

In conclusion, we have presented rigorous results on the diagonalization of a Bogoliubov 
quasiparticle Hamiltonian for a inhomogeneous soliton condensate. Based on the set of 
complete and orthonomalized eigenfunction of the Bogoliubov de Gennes equations, we 
have calculated exactly the quantum depletion of the condensate. We have shown explicitly 
that two degenerate zero-modes, appearing due to the U(l) guage- and a translational- 
symmetry breaking of the system, are the origin of the quantum diffusion and the transverse 
instability of the soliton condensate. We should point out that, in our study presented above, 
a negligible axial trapping potential has been assumed thus the results are valid only for 
a very long cigar-shaped condensate. If the axial trapping potential, which breaks also 
the translational symmetry of the system, can not be neglected and taken as a slowly- 
varying function, a rigorous result on the diagonalization of the Bogoliubov Hamiltonian 
and quantum depletion of a BEC remains to be a challenge. 
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